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In applying the semiempirical intermediate neglect of differential overlap (INDO) method based 
on the Hartree- Fock formalism to a cubic perovskite-based ferroelectric material KNbC>3, it was 
demonstrated that the accuracy of the method is sufficient for adequately describing the small 
energy differences related to the ferroelectric instability. The choice of INDO parameters has been 
done for a system containing Nb. Based on the parametrization proposed, the electronic structure, 
equilibrium ground state structure of the orthorhombic and rhombohedral phases, and F-TO phonon 
frequencies in cubic and rhombohedral phases of KNbOa were calculated and found to be in good 
agreement with the experimental data and with the first-principles calculations available. 



I. INTRODUCTION 

Potassium niobate, a perovskite-type ferroelectric ma- 
terial isostructural to barium titanate, has been subject 
to numerous ab initio electronic structure calculations 
during recent years. Earlier calculations have been per- 
formed for the ideal cubic perovskite structure in or- 
der tOppbtain electron band structure and to interpret 
optical!! or x-ray photoelectrono spectra. Since then, spe- 
cial attention has been paid to the total energy calcu- 
lations making it possible to determine the equilibrium 
geometrytrQ, phonon frequenciesamllj, and interatomic 
interaction-, parameters defining the ferroelectric phase 
transitionsQ'Q 

Most of the calculations cited have been performed us- 
ing the local density approximation (LDA), either with 
the pseudopotential method (Refs. |] and fj]) or with lin- 
earized augmented plane wave (Refs. |, J andjo)) or lin- 
earized muffin-tin orbital (LMTO, Refs. |^ and g)method. 
The latter two approaches use some (different) forms of 
series expansions for the potential inside the muffin-tin 
spheres and in the interstitial. In the LDA-based cal- 
culation schemes, the potential is local and orbital inde- 
pendent, unless Coulomb correlation effects are ad hoc 
introduced by one or .another implementation of self- 
interaction corrections)!^ or within the LDA+U (Ref. 
|l2| ) formalism. This seems to be especially important for 
treating localized states, such as, e.g., those of transition- 
metal impurities in insulators. As has been shown in Ref. 
[l3| for Fe in MgO, the straightforward implementation of 
a LDA scheme may lead to wrong results with respect to 
the energy positioning of impurity levels and the mag- 
netic moment related to the impurity. 

In contrast to LDA, the Hartree-Fock formalism au- 
tomatically incorporates the dependence of the potential 
on the symmetry of a particular orbital, as well as on 
whether this orbital is occupied or not. Another con- 
venient property of the Hartree-Fock formalism that is 
typically realized on a tight-binding basis set is that it 
can be more or less directly applied to crystal surfaces, 
providing correct asymptotics of the electron density at 



the vacuum side. 

Ab initio Hartree-Fock calculations are excessively 
computationally demanding (in the sense that quanti- 
tative results of comparable, or better, accuracy are in 
many cases obtainable within the LDA at much lower 
cost), therefore the applications to perovskite systems 
are not numerous. As an eventually single example, fi- 
nite cluster Hartree-Fock calculation, has been reported 
for a fragment of KNb03 structure. til However, a simpli- 
fication on top of the Hartree-Fock method .known as in- 
termediate neglect of the differential overlapMtJ (INDO) 
lets to decrease the computational effort considerably, at 
the price of treating several parameters as fitting vari- 
ables, to be defined from outside the calculation scheme. 
In contrast to model calculations, which usually require 
an ad hoc fitting, the INDO parameters are believed to 
be largely transferable, so that, once determined for some 
chemical constituent, they may be successfully applied 
in the calculations for a variety of chemical substances 
where the latter participates. 

Typical fields of INDO applications-, include vari- 



s based oru-silicajii 



lie oxides, 



ous defects s\ 

such as— J4ftO Ji§Ej corundum,ML3 zirconiaj^d or alkali 
halides.li3E!l The choice of INDO parameters is not a 
straightforward procedure but rather a trial-and-error 
loop, aimed at reproducing reasonably well band struc- 
ture, equilibrium geometry, and characteristic energy dif- 
ferences for molecules or crystals, as based on experimen- 
tal measurements or ab initio calculations. The list of pa- 
rameters for several elements is given in Refs. |15| and |l| 
along with some discussion on the parameter optimiza- 
tion for ionic crystals. 

The aim of the present paper is to demonstrate that 
the semiempirical INDO method may work well for 
perovskite-type ferroelectrics, and to provide the opti- 
mized set of INDO parameters for all constituents of 
KNb03. Perovskites are generally expected to present a 
problem for any parametrized method, because of varying 
degree of covalency depending on chemical composition 
and because of strong polarizability of transition metal- 
oxygen bonds. An additional difficulty related to fer- 
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roelectric perovskites is that the energy differences that 
play a role in stabilizing the ferroelectric distorted struc- 
ture, due to a fine balance between long-range Coulomb 
forces and short-range chemical bonding, have the or- 
der of magnitude of 1 mRy per formula unit or smaller, 
i.e., much lower than ~1 eV energy differences being d 
cussed, e.g., in relation to charged defects in silica 
Since this is the first, to our knowledge, application of 
the INDO method to perovskite systems, the question is 
to be answered whether the accuracy of the parametrized 
INDO method is sufficient to describe the ferroelectric in- 
stability, and whether the description of the underlying 
energetics is reliable. We optimize the INDO parameter 
set based on the comparison with available ab initio cal- 
culation results and experiments, and answer the above 
question positively by presenting the INDO calculations 
for atomic displacement patterns and phonon frequencies 
that are in good agreement with experimental data. 

The paper is organized as follows. In Sec. 0, we de- 
scribe the essential features of the INDO method and 
the meaning of underlying parameters. In Sec. |l|, the 
choice of the INDO parameters used in our calculation is 
specified, based on the comparison with total energy ab 
initio calculations and the experimental structure data 
for KNbOs. In Sec. IV, atomic coordinates in the 



room-temperature orthorhombic phase and in the low- 
temperature rhombohedral phase are found by the total- 
energy-based structure optimization, and the results of 
INDO calculations for T-TO phonon frequencies are dis- 
cussed. 



II. INDO METHOD AND PARAMETER 
OPTIMIZATION 

The calculation scheme of the Hartree-Fock-Roothaan 
method in the INDO approximation is discussed in de- 
tail in Refs. 15 and |l~6|. Basically, the procedure reduces 
to diagonalizing the matrix of the Fock operator to get 
the one-electron energies, and the linear combination of 
matrix elements with appropriate weights, depending on 
the occupation of corresponding one-electron states, pro- 
vides the total energy. The fixed basis set is minimal 
in the sense that each of the atom-centered functions re- 
lated to the valence-band states (4 in total per oxygen 
atom, 9 per transition-metal atom) is encountered only 
once. The construction of the on-site and off-diagonal 
parts of the Fock matrix is determined in terms of sev- 
eral empirical parameters, labeled by the atom type-|A 
and the index of the atomic orbital (AO) /i (see Rcf.li3). 
The interaction of an electron in the /ith valence AO on 
atom A with its own core 



= -£„cg(M) 



hence Coulomb and exchange integrals and K^, the 
initial guesses for the diagonal elements of the density 
matrix P$} A and for the energy of the /ith AO £^, g (/i), 
i.e., the ion's electronegativity. The interaction of the /ith 
AO with the core of another atom B is approximated as 



V, B 



Z B {I/Rab + - 1/Rab] exp(-aAB-RAs)} , 



contains, apart from the ^ value, which specifies the 
Slater exponent for a one-exponential basis function and 



where Rab is the distance between atoms A and £>, Zb is 
the core charge of atom B, and parameter uab describes 
the non-point character of this interaction. 

Finally, the resonance-integral parameter [3^ enters 
the off-diagonal Fock matrix elements for the spin com- 
ponent u: 

where the /ith and i^th AO are centered at different 
atoms, S^v is the overlap matrix between them, and 
( I ) are two-electron integrals. Parameters (3^ u , a ab 
and E^cgin) are usually fixed throughout the iterations, 

whereas P$ A may be corrected as the self-consistency 
is being achieved. 

It is in principle possible to implement the calculation 
in such way that the diagonalization is done for a num- 
ber of k vectors per iteration. However, conventional 
usage of the INDO method, given the low symmetry and 
possibly the lack of translation invariancc in the systems 
it is usually applied to, restricts the diagonalization to 
the r point of the Brillouin zone only, taking instead a 
supercell, or large unit cell (LUC), all atoms of which 
contribute to the expanded basis set. For ideal systems, 
the enlargement of the unit cell is equivalent to increas- 
ing the density of the k mesh in band-structure calcula- 
tions, since the T point of the reduced (in the supercell) 
Brillouin zone maps onto different points of the original 
Brillouin zone of the single cell. For defect systems, there 
is no problem to treat, e.g., discrete impurity states in 
the dielectric gap, if one or few impurity atoms are in- 
cluded along with tens of bulk atoms in the LUC, since 
the k dispersion of such states is negligible. Anyway, 
the enlargement of the unit cell in case of defect systems 
increases the variational freedom of the basis set. 

Since the construction of the parametrized Fock ma- 
trix is straightforward, the bottleneck of the method in 
what regards the performance and accuracy is the di- 
agonalization of large matrices. Compared to precise 
LDA-based schemes, such as, e.g., full-potential (FP) 
LMTO, which usually pGjwploy multiple-tail representa- 
tion of basis functionsjH3£2l the INDO method manages 
to handle considerably larger supercells. Compared to 
efficient minimal-basis-.computational schemes as, e.g., 
tight-binding LMTO, Id INDO may exhibit such advan- 
tages as unrestricted spatial form of the potential, ab- 
sence of muffin-tin boundary conditions, and of space- 
packing empty spheres. 

As a method essentially based on the Hartree-Fock ap- 
proximation, INDO does not provide a convenient option 
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to incorporate correlation effects into one-electron equa- 
tions, as may be to some extent done within the LDA by 
an appropriate choice of the exchange-correlation poten- 
tial. As a result, the dielectric band gap comes out in 
INDO overestimated usually by 3-5 eV (see, e.g., Ref. |TF 
and 



Moreover, the lack of correlation effects, which 
imply some additional repulsion between electrons, over- 
estimates the chemical binding and results in even more 
underestimated equilibrium bond lengths than is known 
to be the case in LDA calculations. Some part of cor- 
relation corrections (usually referred to as short-range 
corrections as they are mostly of intra-atomic nature) 
may be, however, incorporated in the choice of 1NDO 
parameters, since the latter are based on experimental 
or other external information anyway, and this may to 
some extent improve the two shortcomings mentioned. 



ent properties that are sought to be optimized. 

An example of the total DOS per 135-atom LUC of 
KNbC>3 is shown in Fig. [I] along with the result of FP- 
LMTO calculation. The most obvious discrepancy is in 
the energy separation between the primarily O 2s band 
and the primarily O 2p+Nb Ad valence band. This differ- 
ence is due to the neglect of self-interaction in the LDA- 
based LMTO calculation and the lack of correlation ef- 
fects in the INDO; the experimental x-ray photoelectron. 
measurements set O 2s-0 2p separation at about 15 eV,El 
halfway between the results of Hartree-Fock and LDA 
calculations. The experimental estimate of the optical 
gap of 3.3 eV (Ref. ^5|) is again in between the LDA, 
value of 1.4 eV and 6.1 eV from the INDO calculation.^ 
These differences have a physical foundation and cannot 
be removed without attributing unreasonable values to, 
e.g., O 2s and O 2p-related INDO parameters. 



III. PARAMETER OPTIMIZATION 

In the choice of INDO parameters for our calculation, 
we relied on the experimental information available and 
on the data of ab initio calculations for KNb03 (cited in 
Sec. |), which essentially agree in the description of the 
band structure. Whereas INDO calculations for many 
oxides and potassium salts have been done earlier, and. 
the corresponding parameters for O and K tabulatedi,t£l 
no INDO calculations involving Nb have been, to our 
knowledge, done by now, so one-center and all involved 
two-center parameters had to be found. S ince E ne a is re- 
lated to the central energy position of an AO in question 
which is hybridized with many other states throughout 
the valence band, we calculated partial density of states 
(DOS) by sampling over a single (r) k point in a LUC 
consisting of 2x2x2 or 3x3x3 single perovskite 
cells (40 or 135 atoms in total, correspondingly) and fit- 
ted to corresponding partial DOS from a LMTO calcula- 
tion. 0p,u affects the resonance interaction of the /ith AO 
with other states and hence the width of the correspond- 
ing hybridized bands, which can be as well fitted to the 
a priori known partial DOS. We used for the two-center 
parameter (3^ a weighted value (/3 M + (3 v )/2, therefore 
(3^ and (3 V may be calibrated in such a case as one-center 
parameters. For an initial value of P' ', an expected 
occupation of individual AO's, based on electronegativ- 
ity considerations, in the compound in question may be 
taken, and then refined in the course of iterations. The 
two-center parameter q^b, which does not depend on 
orbital indices, plays a relatively minor role in what re- 
gards the band structure and DOS and affects primarily 
the energetics of atomic displacements, equilibrium bond 
lengths, and hence the equilibrium geometry. Finally, C^, 
nominally being a Slater exponent parameter and as such 
tabulated for all elements, should of course be considered 
here as a free parameter, which is used to improve the 
quality of our fixed, single-exponent basis set. It needs 
some adjustment based on a compromise between differ- 
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FIG. 1. Total density of states of KNbOs calculated with 
full-potential LMTO method (above) and INDO method for 
a 3 x 3 x 3 supercell (below). Energy scale is relative to the 
valence band top. 
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TABLE I. One-center INDO parameters 



Orbital 


C (a.u - 1 ) 


E ncg (eV) 


-P (eV) 


Po (a.u.) 


O 2s 


2.27 


4.5 


16.0 


1.974 


O 2p 


1.86 


-12.6 


16.0 


1.96 


Nb 5s 


2.05 


0.0 


30.0 


0.1 


Nb 5p 


2.05 


-2.0 


30.0 


0.0 


Nb Ad 


1.60 


23.85 


16.0 


0.6 


K 4s 


1.10 


2.8 


2.0 


0.1 


K Ap 


1.25 


0.3 


3.0 


0.03 



The effective charges found by the Mullikan population 
analysis are +0.543 for K, +2.019 for Nb, and -0.854 
for O. This is generally in agreement with simple tight- 
binding calculations ,E3 but emphasizes higher degree of 
covalency of the K-0 bond than may be expected from 
intuitive electronegativity considerations. One should 
note, however, that static effective charges are not well- 
defined properties and vary considerably depending on a 
method used. 

K 3s and K 3p states, which were included into the 
valence-band basis set in the LMTO calculation, have 
been treated as core states in the INDO method. We 
found that in order to obtain correct equilibrium vol- 
ume, it is essential to treat K 3p states as the basis 
AO's within the valence band, since their overlap with 
AO's of other atoms is not negligible. This observation 
is in agreement witluwhat was earlier established in FP- 
LMTO calculations B However, the inclusion of K 3p at 
the expense of K Ap in the minimal one-exponential basis 
of the INDO method does not allow one to describe off- 
center displacements and phonon frequencies with suffi- 
cient accuracy. Therefore, we prefer to keep the K Ap as 
a valence AO and to perform the calculations discussed 
below at the experimental lattice parameters of KNb03. 

Keeping in mind the necessity to obtain reliable val- 
ues of equilibrium atomic displacements and the shape 
of the potential surface related to such displacements for 
subsequent studies of ferroelectric materials, we concen- 
trated on these values as primary criteria for the quality 
of the INDO parametrization we look for. It is known 
that the fine balance between long-range electrostatic 
dipole-dipole interaction and the short-range chemical 
bonding effects is accountable for the ferroelectric insta- 
bility, therefore the parameters cxab and (3^ were espe- 
cially subject to refinement, once E neg and pW are es- 
sentially fixed based on a band-structure analysis. The 
total-energy results from the INDO calculations for dif- 
ferent displacement patterns have been fitted to analo- 
gous data obtained earlier with the FP-LMTO method 
as described in Ref. ||. We made sure that the optimized 
parameter set provides reasonable agreement with the 
FP-LMTO data in describing different displacement pat- 
terns and is not confined to any particular symmetry. 

Since the shape of the total-energy hypersurface over 
atomic displacements is not directly measurable experi- 
mentally, and the results by different ab initio calculation 



schemes differ somehow in determining the depth and 
the position of the off-center potential minima (see, e.g., 
Refs. | and |0|), we relied also on a neutron-diffraction 
data concerning the displaced-atomic positions in the fer- 
roelectric phases of KNb03,c3 and on the V transverse- 
optic (TO) phonon frequencies as additional reference 
points to test our parametrization. The one-center INDO 
parameters we found to provide the best compromise 
in reproducing all these properties are given in Table 
[| The best-fitted two-center parameters cxab are 0.15, 
0.33 and 0.39 a.u." 1 for A=0 and B=0, Nb and K, cor- 
respondingly, and zero for A=Nb and K. The results of 
our ground-state geometry and phonon calculations are 
discussed in the next section. 



IV. RESULTS AND DISCUSSION 

A. Sequence of ferroelectric phases 

As the temperature lowers, KNb03 undergoes a se- 
quence of phase transitions from paraelectric cubic to 
ferroelectric tetragonal then orthorhombic then rhombo- 
hedral phases. The atomic positions in all these phases 
have been determined by Hewat.Ej As a first approxima- 
tion, each of these ferroelectric phases is characterized by 
the off-center displacement of the Nb atom from its sym- 
metric position in the cubic perovskite cell along [100], 
[110], or [111] in three subsequent ferroelectric phases, 
with the gradual lowering of the total energy. On top of 
this major distortion, K and O atoms somehow adjust 
their positions as compatible with the reduced symme- 
try of each particular phase, and a lattice strain even- 
tually appears. The hierarchy of total-energy lowerings 
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FIG. 2. Total energy as a function of off-center Nb displace- 
ments along different directions from its position in the cubic 
perovskite structure as calculated with the INDO (2x2x2 
supercell) method. 
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TABLE II. Calculated T-TO frequencies and eigenvectors in cubic KNbC>3. 



Symmetry 




Eigenvectors (present 


work) 






u> calc 




(cm" 1 ) 






ui es 




(cm" 1 


k 


K 


Nb 





O 


O 


Present 


Ref. 




Ref. 


| Ref. 




Ref. 




Ref. 




Tin 


0.05 


-0.57 


0.70 


0.30 


0.30 


292i 


115i 


203i 143i 


96 




115 


Tin 


-0.88 


0.34 


0.21 


0.16 


0.16 


178 


168 


193 


188 


198 


207 


Tin 


-0.02 


-0.19 


-0.61 


0.54 


0.54 


537 


483 


459 


506 


521 


522 













1 


-1 


272 


266 


234 








280 


a 



a Measurements at 585 K (in the tetragonal phase), Ref. [32] 



related to the Nb displacements along three directions 
is therefore an important benchmark for the quality of 
the calculation in question. In Fig. [^, the energy gain 
due to the Nb displacements from the central position 
in the cubic perovskite cell (with the lattice constant 
a=3.997 A) is shown as calculated by the INDO method 
for the 2x2x2 supercell. 

As is consistent with the experimental data, the [111] 
displacement and hence the rhombohedral phase provides 
the lowest ground-state energy, followed by the [110] dis- 
placement (orthorhombic phase) and the [100] displace- 
ment (tetragonal phase). This qualitative result is rel- 
atively stable against some variations of the INDO pa- 
rameters. As regards the magnitudes of the off-center 
displacements and the depth of the related total-energy 
wells, our INDO parametrization (Table provides good 
agreement with the results of the FP-LMTO calcula- 
tions accounting to all three displacement directions [see 
Fig. 4(a) of Ref. | and Fig. 1(b) of Ref g. 



lated in Ref. g by the FP-LAPW method (for the lattice 
constant a = 4.016 A), and with those calculated in Ref. ^ 
by FP-LMTO. Anyway, the main pattern of atomic vi- 
brations within each particular T\ u mode (primarily Nb 
displacement in the soft mode; almost pure K vibration 
against all other atoms in the intermediate mode, and 
the stretching of the oxygen octahedra in the hard mode) 
are correctly reproduced. The detailed structure of the 
soft mode eigenvector reveals smaller participation of K 
in the displacements with respect to the center of. mass 
than was obtained in the FP-LMTO calculation!] This 
seems to be consistent with the atomic coordinates in the 
tetragonal phase that emerges as the soft mode freezes 
down (see Ref. 28 for the experimental data, and Fig. 1 
of Ref. ||), and this behavior comes out correctly based 
on our INDO parametrization. 



C. Equilibrium displacements in the orthorhombic 
phase 



B. T-TO frozen phonons in the cubic phase 

As another test for the quality of our INDO 
parametrization for the adequate description of the 
atomic-displacement potential surface, we calculated 
the T TO phonon frequencies in the cubic phase of 
KNb03. Similar -calculations have been done earlier 
by other methods ,ot3eI and the experimental dataj-iotj- 
tained mostly by infrared reflectivity measurementsE3H) 
are available. We performed the calculations for a lattice 
constant a = 3.997 A (that is based on an experimental 
perovskite cell volume extrapolated to zero temperature) 
within a conventional frozen-phonon scheme, using the 
2x2x2 LUC. Consistently with the symmetry analysis 
given, e.g., m Ref. §, we studied the effect on the calcu- 
lated total energy of small coupled distortions compatible 
with the T\ u irreducible representation, that reveals three 
TO frequencies, and of the oxygen displacement within 
the single "silent" mode compatible with the Ti u irre- 
ducible representation of the Pm3m space group. The 
calculated phonon frequencies and eigenvectors are given 
in Table ||. 

The calculated frequencies generally fall within the 
limits set by previous ab initio calculations ,t]#0 with 
somehow better agreement for the hard T lu mode and the 
Ti u mode. The eigenvectors agree well with those calcu- 



As an additional benchmark for the fine adjustment 
of two-center INDO parameters, we aimed at obtain- 
ing s~ possibly good agreement with the experimental 
dataEa in determining all atomic positions, and not only 
the Nb displacement, in the orthorhombic and rhom- 
bohedral ferroelectric phases. The orthorhombic phase 
is important because it exists in a broad temperature 
range around room temperature and is subject to most 
studies and practical applications. The rhombohedral 
phase is specially discussed below. Keeping the lattice 
vectors for the orthorhombic phase fixed and equal to 
those listed in Ref. || (a=3.973 A along x = [100] of 
the cubic aristotype, &=5.695 A along y = [011] and 
c=5.721 A along z = [011]), we allowed the c relax- 
ation of K and Nb atoms and the b relaxation of those 
O atoms that are in the same [001] plane with Nb in 
the course of INDO calculations towards self-consistency. 
The total-energy minimization is implemented in the 
code ma king use of the downhill simplex method (see, 
e.g., Ref. |30j ). The resulting atomic p ositi ons within the 
orthorhombic cell are shown in Table III in comparison 



with the neutron-diffraction estimations of Ref. |28|. It 
was of course our aim to provide as good agreement 
as possible by an appropriate choice of INDO parame- 
ters, but the encouraging result is that the agreement 
is very good, given the small number of parameters 
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TABLE III. Positions of atoms in orthorhombic and rhombohedral phases of KNbC>3 (in terms of lattice parameters) as 
determined by neutron diffraction measurements, Ref. Eq, and optimized in the INDO calculation. 



Atom 


a 





c 




A 

^exp 


A 








Orthorhombic phase 








K 








A z 




0.0138±71 


0.0209 


Nb 
Oi 


i 

2 







i 

2 

3 + A Z 




0.0364±10 


0.0347 


On 


1 

I 

2 






A z : 


0.0342±9 


0.0347 


On 


! - A„ 


i + A z f 

4 1 z ) 


A„ : 


-0.0024±9 


-0.0028 








Rhombohedral phase 








K 


A, 


A 2 


A z 




0.0130±81 


0.0139 


Nb 





i 

l+A* 
|+A, 
A z 


i 

§+A* 
A z 

3+ A - 


i 

2 

A, 1 

3 + A, 

|+A, J 


A,, : 
A z : 


0.0301±9 
0.0333±15 


0.0213 
0.0328 



that are not directly related to the structure proper- 
ties. The most noticeable discrepancy is in the relative 
displacement of K atoms, which was somehow overesti- 
mated as compared with the experimental value; the sim- 
ilar trend obtained in the FP-LMTO optimizatioii-of the 
orthorhombic phase was much more pronouncedHil One 
should note that the error in determining the position by 
neutron scattering is, of all atoms involved, maximal for 
K,lj and that our present estimate falls within the error 
bars given in Ref. |2q. 



D. Displacements and phonons in the rhombohedral 

phase 

The rhombohedral phase corresponds to the low- 
temperature ground-state structure of KNbC>3, which 
is what any zero-temperature total-energy minimization 
should normally drive at, with all structural constraints 
lifted. We looked for optimized atomic positions com- 
patible with the symmetry of the rhombohedral phase, 
using the lattice constants of a=6=c=4.016 A (Ref. pg| ) 
but keeping the rhombohedral strain angle fixed. The 
reason for this was that the total energy was found to be 
very insensitive to the rhombohedral strainia BaTiO/3 in 
a FP calculation by Cohen and Krakauerpl and we do 
not expect to achieve better accuracy for KNb03 in our 
INDO calculation. Moreover, we neglected the deviation 
of the rhombohedral strain angle a = 89.83° from 90°. 

The experimental and our optimized atomic positions 
in terms of lattice vectors are given in Table [II. This is, 
to our knowledge, the first optimization of the atomic po- 
sitions in the rhombohedral phase of KNbO/3. The max- 
imal discrepancy with the experiment is for the A x (0) 
parameter that describes a slight stretching of oxygen 
octahedra. This parameter is obviously related to the 
rhombohedral strain and may be slightly adjusted in a 
calculation incorporating the exact value of the strain 
angle. The relative displacements of K, Nb, and O along 



the polar [111] axis are all found to be in very good agree- 
ment with the experiment. The quality of the descrip- 
tion of the total-energy hypersurface in the rhombohe- 
dral phase was further controlled by calculating the T 
TO-phonon frequencies. The general symmetry relations 
between the T phonon modes in the cubic and rhombo- 
hedral phases may be found, e.g., in Ref. |32|. We con- 
sider in the present work only three A\ modes, which 
originate from the T± u block of the cubic phase, as the 
crystal symmetry lowers and all the soft modes become 
stabilized. The c alcu lated frequencies and eigenvectors 
are given in Table IV , The components of the eigenvec- 
tor related to K and Nb displacement exist only along 
[111], whereas each of three equivalent O atoms may also 
have the normal component of the displacement, in the 
threefold axial symmetry along the polar axis. The ex- 
perimental phonon frequency data for the rhombohedral 
phase do not seem to be numerous; the values shown in 
Fig. 8 of Ref. || are about 200, 270 and 600 cm -1 . Our 
calculated frequencies are in good agreement with these 
data. It is interesting to compare the eigenvectors with 
those for the cubic structure. One can see that the "pure 
K" mode is not affected by the structure transformation, 
preserving almost exactly its frequency and the displace- 
ment pattern. The former soft mode of the cubic phase 
only slightly changes the eigenvector, but gets hardened 
up to 278 cm -1 in the rhombohedral structure. Finally, 
the highest-frequency mode has the lowest contribution 
of K and Nb displacements and is essentially related to 
the stretching of the oxygen octahedra, as in the cubic 
phase. 



V. SUMMARY 

In applying the semiempirical INDO method to the 
study of a cubic perovskite system, we demonstrated that 
the method is sufficiently sensitive for the adequate de- 
cription of a ferroelectric instability. The energy gain 
of the order of ~1 mRy per unit cell, i.e., much lower 
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TABLE IV. Calculated frequencies and eigenvectors of the 
T-Ai modes in rhombohedral KNDO3. 



lo (cm x ) 




Eigenvectors 




K||[iii] 


Nb||[iii] 


0||[m] 


Ox[in] 


173 


0.88 


-0.37 


-0.16 


0.04 


278 


0.03 


-0.53 


0.40 


0.28 


593 


0.04 


0.26 


-0.23 


0.51 



than one has to deal with in other conventional applica- 
tions of the INDO method, are nevertheless reliably re- 
produced, resulting in a correct description of the micro- 
scopic structure of ferroelectric orthorhombic and rhom- 
bohedral phases and of the T TO-phonon frequencies and 
eigenvectors. The choice of the INDO parameters was 
proposed for the Nb-containing system, and may be used 
in further applications. 
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